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ABSTRACT 



The theoretical acoustic behavior of an underwater bubble 
screen was examined using Rayleigh reflection theory. A 
microcomputer model simulates the acoustic impedance mismatch 
at the interfaces of an ideal, bubble screen. A sensitivity 
analysis indicates that the angle of incidence of sound 
energy and the speed of sound in the layer are the most 
important screen properties for predicting the insulating 
capability of a bubble screen. In the neighborhood of 
frequencies for which the screen thickness is an integral 
number of half wavelengths, the interference results in a 
reduced reflection coefficient and a corresponding increase 
in transmission through the screen. So that for a broad 
band spectrum wide enough to cover a number of such 
frequencies, the attenuation to be expected exceeds 10 dB 
only over about 90 percent of the spectrum. The interest 
for this work came from the need for a noise insulating 
screen at the Carr Inlet Acoustic Range. 
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I 



INTRODUCTION 



Acoustic trials of ships and submarines are conducted 
Puget Sound Naval Shipyard (PSNS) at the Carr Inlet Acous 
Range in the state of Washington. According to Mro John 
Kriebel of PSNS (personal communication) the measurements 
are often contaminated by noise generated by ship traffic 
in adjacent waters. One method which may offer noise iso 
tion is a bubble screen located between the range and the 
offending noise source. Mr o Kriebel suggested that a 
feasibility study be conductedo 

This thesis is one part of the study and has the 
objective of examining the theoretical benefits of a bubb 
screen, using a computer model based on Rayleigh reflecti 
theory for an homogeneous layer. Recommendations and con 
elusions are made based on a literature search and a 
sensitivity analysis of the computer modelo A companion 
thesis by LT Kelley examines the noise generated by the 
bubbles themselves . 

The writer wishes to acknowledge the interest and 
support from the PSNS Carr Inlet Acoustic Range during 
this study. 
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II 



BACKGROUND 



Aa RESONANCE BUBBLE THEORY 



Experimental studies of the acoustic behavior of bubbles 
in water have been conducted for many years« Normally, when 
a single bubble is being studied, the resonance theory of gas 
bubble pulsations proposed by Minnaert in 1933 |Ref. 1] is 
considered. This theory shows that a gas bubble in water is 
capable of vibrations which become large at an exciting 
frequency, f, in hertz given by 



f = p 



where , 

= ratio of the specific heats of the gas 
P = ambient hydrostatic pressure in pascals 






= density of water medium in kilograms per cubic 
meter 



r = radius of the bubbles in meters 
Table I shows the radius of a resonant bubble as a 
function of frequency f and depth z of concern at Carr Inlet, 



TABLE I 



Resonant Bubble Sizes 



Z 

f 


0 m 




5 0 m 




100 m 


150 


ra 


200 


m 


10 Hz 


0,32 


m 


0 .79 


ra 


1 o 08 m 


1.3 


m 


1,5 


m 


100 


0,032 


m 


0.079 


m 


0 <,108 m 


0.13 


m 


0.15 


m 


1000 


3 . 2 


ram 


7,9 


mm 


1 0 <, 8 ram 


13 


mm 


15 


mm 


10,000 


0.32 


mm 


0.79 


mm 


1 0 8 ram 


1,3 


mm 


1.5 


mm 




wher 




‘ = 1-? 


)3 X 


10^ kg/ra 


3 












P 


= 10" 


' Cl 


+ 0<,lz) 


pascal 









^ = 1.40 
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Carstenson and Foldy [Refo 2] examined certain acoustic 
properties of bubble screens and demonstrated reasonable 



experimental agreement with the 
bubbles in water. 

Although resonance theory p 
tions of incident sound energy 
attenuation is a very strong fu 
"75 percent of the energy lost 
bubbles of radii within 10 perc 
[Refo 3: pol3j A parameter us 

is its extinction cross section 
energy lost from the sound wave 
the bubble. [Ref. 4j ^ is mad 
section (OT ) and scattering cro 



theory of resonance of 

redicts very large attenua- 
for bubbles in water, this 
notion of the bubble sizeo 
by the sound wave is due to 
ent of the resonant radiio" 
ed to describe bubble behavior 
(^) , the ratio of the 
to the intensity incident on 
e up of both absorption cross 
ss section (<JZ ) where 

o 



eg: = <J7 

For a plane wave propagating in a bubbly mixture, the 
attenuation, a, for bubbles of uniform radius r is 
{Ref . 5 : p . 6 15J 

a = Ik dB per meter 

where n(r) = number of bubbles per cubic meter of 
radius r 

Medwin [Ref« 6j has indicated that "^oc,at a resonance, the 
scattering and absorption cross sections of a typical bubble 
at sea are of the order 1000 times its geometric cross 
section o '* 
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Several points should be made here. 

1. The sensitivity of attenuation on bubble sizes 
requires a precise control of bubble size. 

2. As will be discussed later, even if bubbles of 
radius greater than 1 to 2 cm could be produced, these 
bubbles will rapidly become unstable and break up. 

3. A bubble will expand mono tonically as it rises to 
accomodate the decrease in hydrostatic pressure. A bubble 
of 1 mm radius at 100 m will expand to 3.2 mm at the surface 
and its resonance frequency will vary from 10 kHz to 1.0 

kHz during its ascent. 

4. MacPherson [Ref. 7J indicates that it is difficult 
to produce bubbles of a given size. Consequently, it is 
difficult to predict the attenuation obtainable from given 
bubble production mechanisms. 

B. RAYLEIGH REFLECTION THEORY 

For the above reasons and because of the difficulty of 
producing bubble resonance at low frequencies, this paper 
proposes to view the bubble screen as an homogeneous layer 
separating semi-infinite layers of bubble free water. 
Rayleigh reflection theory predicts that attenuation of 
sound energy by the layer of bubbles will be a function of 
layer width, angle of incidence, frequency, and the acoustic 
properties of the bubble layer. A layer of constant width 
will be used to simulate the bubble screen. 
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A sensitivity analysis will then indicate which 
parameters are of greatest concern in bubble layer design. 

The sound power transmission coe f f i c i en t , , will be the 

measure used to analyze these parameters. is related to 

the sound power reflection coefficient 0\ by the equation, 

= 1 

C. MANUFACTURING BUBBLES 

In the past, researchers used a variety of techniques to 
control the size of bubbles in screens. The methods attempted 
include electrolysis, J_Ref. 8j microdispersers, [Ref. 9j 
porous porcelain filters, [Ref. lOj and perforated pipes. 

Exact bubble size was important to determine acoustic 
characteristics of scattering cross section, extinction cross 
section, resonance frequency, and other parameters. Several 
papers concerning bubble generation [Refs. 11, 12, 13, 14] 

indicate that some of the many factors controlling bubble 
size for perforated pipes include: 

1 . orifice size 

2. pressure drop across the orifice 

3. wetting of the orifice 

4. contaminants in the media 

5. angle of inclination of the orifice 

6. rate of bubble formation 

Uncertainties in these factors are the primary factors 
which will prevent accurate bubble size prediction for a 
bubble screen installation on a large scale such as at Carr Inlet. 
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BUBBLE HYDRODYNAMICS 



Air bubbles in water exhibit complex hydrodynamic 
characteristics. However, their behavior can be categorized 
generally with respect to the bubble size. Initial bubble 
motion has been analyzed by Walters and Davidson 
[Refs. 15 and 16j and Bachhuber and Sanford jRef. 17j. 

Gas diffusion from bubble affects bubble behavior signif- 



icantly 


when th 


e 


bubble diameter 


is small. 


The 


gas diffusion 


effect 


depends 


on 


the rate of asc 


ent, the 


radius 


of the 


bubble 


and the 


di 


fference between 


the gas 


concen 


tr at ions 



across the bubble liquid interface. [Ref. 18: p. 714j 



There are 


two critical radii: 








R 


the radius above which 


the 


bubble 


grows mono to ni cal ly 


a 


in size as it ascends; 


gas 


leakage 


is not predominant 




the radius below which 


the 


bubbles 


shrink and 




collapse bubbles whose 
Rb are unpredictable. 


rad 


i i lie b 


etween Ra and 



Two bubbles of the same gas and size produced at the same 
depth can be affected differently by surfactants, causing 
one of the bubbles to grow and the other to shrink and 
collapse. 

Bubbles of interest for a screen at Carr Inlet are in 
the 1-10 mm diameter range. This size has proven to be the 
easiest to produce via a perforated pipe. This method is 
also by far the most convenient and cost effective. 

"Hydrostatic head has very little influence on the 
relation between the rate of formation and the size of the 
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bubble . ** 



[Ref. 19: p.l7j VanKrevelen and Hoftijzer 
[Ref. 20: p.30j showed that "...the diameter of the bubbles 

is independent of flow rate and increases with the cube root 
of the orifice diameter." This statement breaks down at 
higher flow rates when gas jetting effects begin. 

Crump [Ref. 2l] suggests a differential pressure of 2-3 
psi be established between pipe line pressure and hydrostatic 
water pressure. This should make the desired bubbles for 
forming the bubble screen. A bubble of radius Imra will 
expand mono toni cal ly due to decreased hydrostatic pressure 
and will not be affected to a significant amount by gas 
diffusion and temperature effects. A 1 mm diameter bubble 
at 150 m depth will increase in size to 2.5 mm in ascending 

isothermally to the surface. 

A sample calculation of this effect is shown in Table 

II and shows that this bubble will shift its resonant 

frequency by several kilohertz during its ascent. 

TABLE II 

Typical Bubble Size Expansion 

RADIUS in mm RADIUS in mm RESONANCE FREQUENCY in kHz 

(600 f t . ) ( surface) (600 ft.) (surface) 



. 200 


. 5 


70.4 


16.0 


. 500 


1.23 


28.2 


6.4 


1.0 


2.5 


14.1 


3 . 2 


2 . 0 


4.9 


7 . 0 


1.6 


10.0 


24.6 


1.4 


.13 
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RISE OF BUBBLES 



E . 



The rate of rise of bubbles has been thoroughly 
documented. jRefs. 22, 23, 24, 25] The speed at which a 

bubble rises depends primarily upon the size and shape of 
the bubble, and can be categorized as in Table 



[Ref. 26: p. 72] 



TABLE III 

Characteristics of Rising Bubbles 
SHAPE DIAMETER in MM SPEED in cm/s Re NUMBER 



MOTION 



spherical 



.5 



5-10 



< 200 rectilinear 



ellipsoidal 1-15 



spherical- >15 

capped 



10-30 200<Re<4700 zig-zag 

helical 

oscillating 

30-40 >4700 erratic 

unstable 



Increasing bubble size causes: 

1. increased coefficient of drag 

2. increased vertical speed 

3. nonlinear motion 

4. oscillating motion 

5. wake shedding 

Turbulence seems to be the controlling factor in the motion 
of bubbles larger than 1,200 mm. [Ref. 27: p. 121] "Bubbles 
smaller than about 0.1 cm in water rise along smooth paths 
and do not generate any measureable sound when rising." 

[Ref. 28: p. 25] 
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Surface contaminants affect bubble behavior in many ways. 
These surfactants [Ref. 29] will have their greatest effect 
on the terminal speed of a bubble by affecting transition 
from rectilinear motion to oscillating motion. The bubble’s 
wake structure is changed causing delayed boundary layer 
separation, delayed vortex shedding and smaller bubble wakes. 

Terminal speed will be a function of bubble diameter, 
surfactants and Reynold’s number. Approximate increase in 
vertical speed due to multiple bubbles range up to a factor 
of 1.5 times the single bubble velocities. [Ref. 30] 

F. BUBBLE NOISE 

Previous investigations [Ref. 31] show that greater than 
90 percent of sound produced by bubbles in water occurs during 
bubble formation. 



o o n o o 




Figure 1. Sound produced by a bubble at an orifice 



This sketch [Ref. 32: p. 23] correlates the process of an 
individual bubble leaving an orifice with the observed time 
variation of the radiant noise. The peak sound pressure 
occurs at the time the bubble separates from the orifice. 

Strasberg [Ref. 33: p.24] predicted the peak sound 
pressure at distance d to be. 
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I?n,s 




V2 



total dissipation constant 
•014 for an air bubble in water 
total volume rate of bubble formation 
excess pressure 
This works out to be, 

- 132 dS re 1 ^Fa at 1 m. 

for a bubble rate of 10/sec for f = 1.0 kHz 
After separation, larger bubbles generate continuous 
sound levels of relatively small amplitude {46 to 66 dB re 
1 /iPa at Im). [Ref. 34: p.25] 

This is apparently due to oscillations induced by their 
irregular flow path. Smaller bubbles produce less or 
immeasurable sound when rising. The coalescing of bubbles 
cause one tenth the sound of bubble formation. 



where , 



cl = 

V = 

p = 



G. PROPERTIES OF BUBBLE SCREENS 

Studies on the acoustic theory of bubbles and bubble 
screens were done during World War II. [Refs. 35, 36, 37] 

Efforts in the United States were concentrated through the 
National Defense Research Committee (NDRC) under a broad 
research program supported by the Office of Scientific 
Research and Development (OSRD) . The prominent NDRC 
organizations whose findings were published by the U.S. 
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Navy after the war as NDRC Summary Technical Reports were 
the University of California Division of War Research 
(UCDWR ) , the Columbia University Division of War Research, 
the Columbia University Sonar Analysis Group, and the Woods 
Hole Oceanographic Institution. The UCDWR was the largest 
of the associated groups and came to be known as the Sonar 
Data Division, accounting for a significant amount of 
empirical studies on ship's wakes, scattering and absorption 
of bubbles, and sound transmission through bubble screens. 
Carstenson and Foldy were part of this effort and in 1947 
they published their study of the propagation of sound 
through a bubble screen. Their screen was approximately 
17 in. long and 10 ft. high with thickness varying from 
3 in. to 6 in. The pulsed screen was analyzed as if 
bubbles were uniformly distributed between two infinite 
parallel planes and the speed of sound in the screen was 
assumed to be approximately the same as that of water. 

For the case of continuous - flow screens, the bubbles 
were of varying sizes and speed of sound in the layer 
considerably different from that in water. 

Their results supported the resonant bubble theory for 
pulse screens at frequencies from 15 kHz to 35 kHz. However, 
their data from 5 kHz to 45 kHz for continuous flow screens 
of bubbles of varying sizes were not very satisfactory. 
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Throughout their study, the thickness of the screen was 
assumed to be the arithmetic average of the observed layer 
width. The density of the screen was based on observations 
of the average number of bubbles present per unit volume 
and of the average bubble size. 

For the continuous flow case, the screen had a "..fairly 
sharply defined core of larger bubbles ... but in front and 
back there is a gradually tapering distribution of smaller 
bubbles." [Ref. 38: p.SOlJ As a result, the speed of sound 

was a gradually changing function of the distance normal to 
the layer interface plane causing reduced reflections of 
sound . 

Wood and Spitzer [Refs. 39,40j proposed theoretical sound 

speeds which would occur for given concentrations of air by 

volume in water. Figure 2 illustrates those results and 

experimental verification by others. [Refs. 41, 42] For 

-2 -3 

concentrations between 10 and 10 , readily obtainable in 

the laboratory or in the field, sound speeds as low as 
100 - 300 m/sec are found. 

H. APPLICATIONS OF LAYERS 

There are extensive studies of the reflection and trans- 
mission of energy in layered media. [Refs. 45, 46] Numerous 

applications of the use of layers range from reflection 
reduction for optical lens coatings to absorbent materials 
for architectural acoustics. The ratio of layer thickness 
to the wavelength of the energy in the layer is a very sig- 
nificant parameter due to the standing wave pattern produced. 
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Figure 2. Sound speed in air-water mixtures 
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III. BUBBLE LAYER AT CARR INLET 



A. OCEANOGRAPHIC CHARACTERISTICS 

The chart of Carr Inlet Acoustic Range in Appendix A 
indicates a large, irregularly shaped area surrounding the 
test site. The site is roughly 5 nm long and 3 nm wide. 

It is proposed to consider in this study a bubble screen 
crossing the mouth of the inlet, a distance of approximately 
2900 m. This position places the screen directly between 
commercial boat traffic to the southeast and the test site. 

The water depth at this position varies from approximately 
450 ft. at the relatively level bottom at the middle of the 
inlet, to relatively steep shore gradients approaching Fox 
Island and McNeil Island as shown in Figure 3. 

The frequencies of interest for attenuation will be 
broadband (100 Hz to 10 kHz+) . 

Oceanographic studies [Ref. 47j at Carr Inlet indicate 
that the tides and currents which are likely to be encountered 
will be variable but less than one-half kt. There is no 
significant evidence for underground springs or bottom 
currents which would cause any significant salinity or 
temperature gradients. Seasonal variations in water 
temperature and salinity occur due to the annual air tem- 
perature cycles and rain runoff. However, no significant 



CARR INLET - SIDE VIEW 



LCXX> 



'^Trrrr ^ r// rrr 7h-7'rTT--r7 rj > V7///7 /// 



f screen 



1 in = 1 nm 



depth 
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Figure 3. Position of bubble screen in sensitivity analysis. 
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impact is anticipated on the sound velocity profile. So, 
for the distances encountered in this application, the 
sound speed in water will be assumed to be constant 
(C = 1500 m/s ) . 

B. SCREEN POSITION 

Considering these oceanographic factors, the bubble 
screen will be assumed to be produced at the bottom of 
the inlet and to rise vertically with constant width. An 
estimation for actual layer width spreading will be 
presented in Chapter V. 
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IV. SENSITIVITY ANALYSIS 



A. THEORY 

At frequencies removed from the resonant frequency of 
the bubbles, a bubble layer is expected to reflect sound as 
a result of two interrelated effects on the acoustic 
properties of the medium. [Ref. 48: P.925] First, the 
compressibility of the water is changed with, the introduction 
of the bubbles which in turn changes the sound speed in the 
layer. Second, density will be changed. For the bubble con- 
centrations considered in this thesis, this latter change is 
small and the principal acoustic characteristic of the bubble 
layer will be its sound speed. The sound speed in bubble 
layers was addressed by Wood [Ref. 49j whose work was con- 
firmed experimentally by Silberman. [Ref. 50: p.90] 

In the present work, the acoustic effect of a bubble 
screen will be analyzed in terms of the sound power trans- 
mission coefficient using a computer simulation of the bubble 
screen. As indicated earlier, this computer simulation will 
be based on Rayleigh reflection theory for an homogeneous 
layer of bubbly water with constant thickness. 

Rayleigh reflection theory is based on a solution of a 
system of acoustic equations, describing the conditions that 
must be satisfied at each boundary of the layer. 
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The boundary conditions at the interfaces are continuity 
of the acoustic pressure and of the normal component of the 
particle velocity. 



For oblique incidence and when the media in front and 
behind the layer are the same, the sound power transmission 
coefficient C><^ reduces to [Ref. .51], 







This coefficient represents the ratio of the intensity of 
the transmitted sound to the intensity of the incident sound. 
The geometry of the problem is shown in Figure 4, and 
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Figure 4. Bubble layer reflection geometry 
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AX 



is a function of width and acoustic impedance of the 
layer and depends on the angle of incidence and frequency of 
the incident sound. 

The subsequent sensitivity analysis will indicate the 
theoretical performance expected when parameters vary. 

B. SENSITIVITY ANALYSIS RESULTS 
1 . Parameter description 

Data were taken using the options in Table IV. A 
change in symbols was required for use in the computer model. 
The new symbols used to represent the parameters are also 
listed for the options in Table IV. 

TABLE IV 

Analysis Options 

OPTION ORDINATE ABSCISSA CONSTANT PARAMETER 
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Each of these options will be addressed in the following 
analysis. Values of A and their corresponding fraction 
volume of air bubbles are shown in Table V, 

TABLE V 



Sound speed ratios for given air concentrations 



A 


Fraction 

air 


by 

in 


volume of 
water 


in 
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8 


X 


io“^ 


0 . 2 


2 


X 
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These values 
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be 


expected in 



practical bubble screen. 

The value of L2 was obtained by incorporating a normal- 
ization factor which forced the data to repeat itself after 
L2 reached a value of 1.0. This was done for ease in 
plotting and clarity in data presentation. The normalization 
factor LI is found by noting that in Eq. (1) , 



OK' — 1 ^ for Id^ a — ^ TTT n — Oj 1 ^ 2, ^ , 



for what values of layer width/wave length ^ j? 






iLlh 



cos 



a = n7T 



C- 



31 



Let n 



1 



whe 

The 

2 . 

run the 
"0PT6" , 
tion of 
program 



- h 



-^ 4 “ ^ “ 

; 



re LI = and, by Snell’s 



re 



fore, LI = 1 1 - 5 



The Computer Program 
The WANG SYSTEM 2200S desk 
analysis and provide data o 
was written in BASIC langua 
the variables and svmbols a 



The plotting option is sele 
executed. Each option has its own 
interrogates the operator for the b 
required to produce the selected da 
for all the parameters are pre-sele 
transmission coefficient is plotted 
changing option parameters. 

3 . Option 1 

When the sound power transm 
plotted against the normalized laye 
tween maxima and minima. Minimum s 
through the layer at odd integermul 
wavelengths. Figures 5 to 10 show 
L2 of a fixed Al for a variety of v 



1 

Law, cc, 

and L2 = — 

^< 2 . 

top computer was used to 
utput. The program, 
ge and provides a defini- 
t the beginning of the 

cted when the program is 
subroutine which 
ubble screen parameters 
ta output. Step sizes 
cted. The sound power 
vertically against the 

ission coefficient is 
r ratio, T3 cycles be- 
ound power is transmitted 
tiples of one-quarter 
the dependence of T3 on 
alues of A . 




32 




A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 

Figure 5. Sound power transmission coefficient vs. 

layer ratio, parameter A for constant 
Al = 0° 
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A = SOUND SPtED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 

Figure 6. Sound power transmission coefficient vs. 

layer ratio, parameter A for constant 
A1 = 20° 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 7 , 



Sound power , 
layer ratio, 
A1 = 40° 



transmission coefficient vs. 
parameter A for constant 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 8. Sound power transmission coefficient vs. 

layer ratio, parameter A for constant 
A1 = 60° 
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A = SOUND SPttiD RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ wavelength RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figur 0 9. Sound. powGr# transmission cosfficisnt vs, 
layer ratio ^ parameter A for constant 
A1 = 80° 



) 
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T3 IjOr 




A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH /WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 10. Sound power transmission coefficient vs. 

layer ratio, parameter A for constant 
Al = 85° 
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The following observations can be made: 



a. Regardless of angle of incidence, if the sound 
speed ratio is less than 0.1, attenuation of at least lOdB 
can be expected over at least the middle 75 percent of the 
layer ratio range. 

b. Regardless of the sound ratio or the angle of 
incidence, there will always be at least 5-10 percent of the 
layer ratio which will have only negligible attenuation. 

4 . Option 2 

When the sound power transmission coefficient is 
plotted against the angle of incidence, the maximum trans- 
mission always occurs at normal incidence and the minimum 

o 

transmission at angles approaching 90 . Figures 11 to 16 

are plots with the layer ratio held constant and the sound 
speed ratio varied as a parameter. 

The following observations can be made: 

a. When the sound speed ratio reaches 0.1, at 
least lOdB attenuation can be expected over all angles of 
incidence and over at least 75 percent of the layer ratio 
range . 

b. The only values of sound speed which will give 
a uniform value of attenuation over all angles of incidence 
are speed ratios of less than 0.1. This situation would 
require the bubble screen to maintain a volume fraction 

of air to water of greater than approximately 10 
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T3 0 



30 " 60 « 90 ” 

A1 

A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 11, Sound power transmission coefficient vs, 
incident angle^ parameter A for constant 
L2 = 0 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 12. Sound power transmission coefficient vs. 

incident angle, parameter A for constant 
L2 = . 1 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 13. Sound power transmission coefficient vs. 

incident angle, parameter A for constant 
L2 = . 2 
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T3 




A1 

A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 14. Sound power transmission coefficient vs. 

incident angle, parameter A for constant 
L2 = . 3 
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A1 



A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 15. Sound power transmission coefficient vs. 

incident angle, parameter A for constant 
L2 = . 4 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 16. Sound power transmission coefficient vs. 

incident angle, parameter A for constant 
L2 = .5 
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5 . Option 3 



As in the plots for option 2 , the sound power 
transmission coefficient is plotted against the angle of 
incidence. This time, though, the sound speed ratio is 
plotted as a parameter. See Figures 17 to 22. 

The maximum sound transmitted occurs at normal 
incidence and the minimum always occurs as the angle of 
incidence approaches 90^. 

The following observations can be made: 

a. The significance of a change in the sound 
speed ratio is very clear in this option. If Figure 18 
(sound speed ratio = 0.1) is compared to Figure 22 (sound 
speed ratio = 0.5) , a dramatic decrease of attenuation is 
seen to occur. 

b. At low angles of incidence, values expected for 
attenuation are relatively constant. At high angles of 
incidence, attenuation values are changing rapidly, 
indicating less confidence for a given predicted layer width. 

6 . Option 4 

This option plots the sound power transmission 
coefficient against the normalized layer ratio. The angle 
of incidence is the parameter which is varied for constant 
values of sound speed ratios (Figures 23 to 29) . Again it 
is seen that values for the sound power transmission 
coefficient cycle between maxima and minima. The minimum 
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T3 ^ 



L2 = .i 




A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 17. Sound power transmission coefficient vs. 

incident angle^ parameter L2 for constant 
A = . 05 
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A = SOUND SP(£ED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 18. Sound power transmission coefficient vs. 

incident angle, parameter L2 for constant 
A = . 1 
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T 3 ’-0 




A = SXDUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 19. Sound power transmission coefficient vs. 

incident angle, parameter L2 for constant 
A = .2 
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T 3 i.of 




A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 20. Sound power transmission coefficient vs. 

incident angle, parameter L2 for constant 
A = . 3 
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A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 21. Sound power transmission coefficient vs. 

incident angle, parameter L2 for constant 
A = . 4 
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T3 1.0 




A1 



A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 22. Sound power transmission coefficient vs. 

incident angle, parameter L2 for constant 
A = .5 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 23. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 
A = .05 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 24. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 
A = .1 
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A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 25. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 
A = . 2 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



pj_gu]T6 26. Sound power transmission coefficient vs. 

layer ratio, parameter Al for constant 
A = . 3 
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A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 27. Sound power transmission coefficient vs. 

layer ratio, parameter A1 for constant 
A = . 4 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 28. 



Sound 
layer 
A = . 



power 

ratio 
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transmission coefficient vs. 
parameter Al for constant 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 29. 



Sound 
layer 
A = . 9 



povrer 

ratio 



transmission coefficient vs. 
parameter Al for constant 
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values occur at layer ratios which correspond to odd 
multiples of one-quarter wavelength. 

The following observations can be made: 

a. Changes in attenuation due to variations in 
the angle of incidence are very slight for angles less 
than approximately 45^, but the attenuation in uniformly 
predictable and by far the greatest at angles greater 
than 80°-85°. 

b. Variations in the sound speed ratio cause far 
greater effects on the attenuation when the sound speed 
ratio is greater than 0.2. 

7 . Option 5 

The sound power transmission coefficient is plotted 
against the sound speed ratio in this option. The bubble 
layer ratio is kept constant while the angle of incidence 
is varied. A much different view of the data is seen in 
Figures 30 - 34. 

The following observations can be made: 

a. The best attenuation occurs at very low sound 
speed ratios as was also seen in options 1 and 2, 

b. There is a significant amount of attenuation 
gained for angles of incidence past 75^. 

8 , Option 6 

In this option, the sound power transmission 
coefficient is plotted against the sound speed ratio with 
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the parameter of layer ratio varied for constant values 
of angles of incidence in Figures 35 to 38. 

The following observation can be made: 
a. Layer width ratio curves vary only slightly 
indicating that layer width is not an important parameter. 
In addition, there is still approximately 10 percent of 
the layer ratio range (L2 < .05) virtually unaffected by 
sound speed ratio or angle of incidence. 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 30 . 



Sound power 
sound speed 
L2 = .05 



transmission coefficient vs. 
ratio, parameter A1 for constant 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 31. Sound power transmission coefficient vs. 

sound speed ratio, parameter Al for constant 
L2 = .1 
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•5 1.0 

A 



A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH / WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 32. 



Sound power 
sound speed 
L2 = . 2 



transmission coefficient vs. 
ratio, parameter Al for constant 



6ft 




A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 33. Sound power transmission coefficient vs. 

sound speed ratio^ parameter Al for constant 
L2 = . 3 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 34, Sound power transmission coefficient vs. 

sound speed ratio^ parameter A1 for constant 
L2 = , 4 



) 



66 




A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 35. Sound power transmission coefficient vs. 

sound speed ratio, parameter L2 for constant 
Al = 0° 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 36. 
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transmission coefficient vs. 
ratio, parameter L2 for constant 
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T3 1-0 




A 

A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 



Figure 37. Sound power transmission coefficient vs. 

sound speed ratio, parameter L2 for constant 
A1 = 60° 
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A = SOUND SPEED RATIO 

A1 = INCIDENT ANGLE 

L2 = LAYER WIDTH/ WAVELENGTH RATIO 

T3 = SOUND POWER TRANSMISSION COEFFICIENT 
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transmission coefficient vs* 
ratio, parameter L2 for constant 



Figure 38, 
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BUBBLE SCREEN APPROXIMATIONS AND RECOMMENDATIONS 



A. FREQUENCY BANDS WITH LITTLE ATTENUATION 

Options 1 and 6 shows that there will be at least 10 
percent of the layer ratio range virtually unaffected by 
a bubble screen. So for any given spectrum, there will be 
frequency bands of small attenuation which occur when the 
screen thickness is an integer multiple of one-half 
wavelengths in the layer. 

B. LAYER SPREADING 

Answers to questions on layer-width spreading were not 
found in the literature. Forcing parameters which are 
present in Carr Inlet which would cause a bubble screen to 
lose its simulated parallel interfaces with the surrounding 
media are bubble interactions, turbulent bubble motion, 
and inlet currents. In order to understand these effects, 
it might be better to conduct an experimental study at 
Carr Inlet and not a theoretical study because predictability 
rapidly disappears due to turbulent effects when bubbles 
grow larger than 1 to 2 cm. 

However, an estimate of the layer spreading may be 
taken from a linear extrapolation of the observed laboratory 
results of Carstenson and Foldy. jRef. 52j Their screen 
spread up to 6 in. for a rise of 10 ft. Figure 39 shows 
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Figure 39. Bubble screen spreading effect 
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RECOMMENDATIONS AND ALTERNATIVES 



D . 

The purpose of this thesis w^as to investigate the 
benefits of a bubble screen which may offer some isolation 
from underwater noise interference at Carr Inlet, As a 
result of this sensitivity analysis and further investiga- 
tion^ the following recommendations are made: 

1. An alternate screen position as shown in Figure 40. 




Figure 40. 



Alternate bubble screen position 



a. The advantages of this position are; 

CD Only a fraction of the piping (one estimate 
for piping to cross the inlet for a uniform screen required 
either a large number of in line pressure regulators or 
almost 100 km of piping) . 

(2) Shallower depth allowing less dispersion of 
bubble layer. 

(3) Exploitation of the high attenuation pre- 
dicted in this sensitivity at steep angles of incidence. 

b. The disadvantages include: 

(1) Increased noise interference when screen 
is placed closer to the test site. 

(2) Low frequency sound still not affected. 

C3) Diffraction of sound around edges reduces 

any attenuation. Further study for any perforated pipe 
bubble screen mechanism should include pipe suspension or 
anchoring devices, pip^ cleaning methods and compressor 
r equir erne nt s . 

2. Conduct an on-site layer spreading experiment. 

3. As an alternative to the bubble screen solution, 
investigate the benefits of adaptive beam forming to null 
out bothersome interference directions. 
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VI. CONCLUSIONS 



A. The insulating capability of a bubble screen at Carr 
Inlet is a function of the angle and frequency of incident 
sound energy and the width and speed of sound of the bubble 
layer. Of these parameters, the angle of incidence and the 
speed of sound in the layer are the most important to consider. 

B. There will be approximately 10 percent of the frequency 
spectrum over which virtually no attenuation of sound will 
occur no matter what speed of sound could reasonably be 
acheived in the bubble layer. 

C. The maximum consistent attenuation which could be 
expected is approximately 10 dB at bubble concentrations 
greater than about 5 x 10 

These conclusions are based on the assumption that the 
air-water mixture in the layer is an ideal, homogeneous fluid. 
As a next step in the improvement of this model, the case 
of a lossy medium should be incorporated to include the 
effects of absorption of a bubble screen. 
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APPENDIX A: BUBBLE SCREEN POSITION AT CARR INLET. 
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COMPUTER PROGRAM 



10 REM ; PROGRAM NAME "OPT6" 

20 REM ; PURPOSE ... PLOTS THE POWER TRANSMISSION COEFFICIENT 
(T3) FOR A PLANE WAVE INCIDENT ON A BUBBLE SCREEN LAYER IN 
WATER. THE DESIRED OUTPUT DISPLAY IS SELECTED FROM ONE OF 
THE OPTIONS BELOW 

30 REM ; DEFINITION OF VARIABLES 

40 REM ; "A" IS ACOUSTIC IMPEDANCE RATIO = R2/R1 
50 REM ; "Al" IS THE INCIDENT ANGLE MEASURED FROM NORMAL 
60 REM ; "A2" IS THE ANGLE OF TRANSMISSION IN THE SCREEN 
70 REM ; "C" IS THE SOUND SPEED RATIO=C2 / 15 00 = "A" APPROX 
80 REM ; "C2" IS THE SOUND SPEED IN THE BUBBLE LAYER 
90 REM ; "L" IS THE RATIO, BUBBLE SCREEN WIDTH/WAVELENGTH 
100 REM ; "LI" IS THE NORMALIZING FACTOR FOR L 
110 REM ; "L2" = L/Ll 

120 REM ; "T3" IS THE SOUND POWER TRANSMISSION COEFF (OFFICER ) 
130 REM ; "S" IS THE SELECTED PLOTTING OPTION 

140 REM ; PLOTTING OPTIONS 

150 REM ; OPTION 1 IS T3 VS L2 FOR CONST Al AND PARAMETER A 

160 REM ; OPTION 2 IS T3 VS Al FOR CONST L AND PARAMETER A 

170 REM ; OPTION 3 IS T3 VS Al FOR CONST A AND PARAMETER L2 

180 REM ; OPTION 4 IS T3 VS L2 FOR CONST A AND PARAMETER Al 

190 REM ; OPTION 5 IS T3 VS A FOR CONST L AND PARAMETER Al 

200 REM OPTION 6 IS T3 VS A FOR CONST Al AND PARAMETER L2 

210 INPUT "OPTION NUMBER =",S 

215 SELECT R 

220 IF S = 1 THEN 300 

230 IF S = 2 THEN 570 

240 IF S = 3 THEN 840 

250 IF S = 4 THEN 1100 

260 IF S = 5 THEN 1380 

270 IF S = 6 THEN 1700 

28 0 REM ; SUBROUTINE OPTION 1 

290 REM ; THIS SUBROUTINE PLOTS T3 VS L2 (Al CONST), PAR A 
300 INPUT "INCIDENT ANGLE IN RADIANS(EA 10 DEG=.174 RAD)=",A1 
310 SELECT PRINT 213(64) 

320 PRINT " INCIDENT ANGLE=",A1 

330 PRINT " T3(0-1.0) VS L2(0-1.0), PARAMETER A(.I-.9)" 

340 SELECT PRINT 413 

350 PLOT (10) (0,-50,"+") 

360 PLOT (10) (50, 0, "+") 

370 PLOT (-500, 0,) 

380 FOR A = .1 TO .9 STEP .4 
390 H=0: V=0: H1=0: V1=0 
400 A2 = ARCSIN(A*SIN(A1) ) 

410 Ml = A"COS(Al)/COS (A2 ) 

420 M2 = (l/A)*COS(A2)/COS(Al) 

430 LI = 1/(2*SQR(1-A*A*SIN(A1)*SIN(A1) ) ) 

440 FOR L2 = 0 TO 1 STEP .02 
450 B2 = 2*#PI*L2*L1*C0S(A2 ) 
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460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 

610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 



D1 = COS(B2)‘‘COS(B2) 

D2 = SIN(B2)*SIN(B2) 

T3 = l/(Dl+( (M1+M2)*(M1+M2)/4*D2) ) 

H = INT(L2*500): V = INT(T3*500) 

PLOT( (H-Hl) , (V-Vl) ,HEX(FB) ) 

HI = H: VI = V 

NEXT L2 

PLOT ) 

NEXT A 
GOTO 210 

REM ; SUBROUTINE OPTION 2 

REM ; THIS SUBROUTINE PLOTS T3 VS A1 (L CONST), PARAM . A 
INPUT "BUBBLE LAYER/WAVELENGTH= " , L 
SELECT PRINT 213(64) 

PRINT " BUBBLE LAYER/WAVELENGTH= " , L 

PRINT " T3 (0-1.0) VS AKO-89 DEG), PARAM A(.l-.9) 

SELECT PRINT 413 

PLOT (10) (0,-50,"+") 

PLOT (10) (50,0,"+") 

PLOT (-500,0,) 

FOR A = .1 to .9 STEP .4 
H=0: V=0: H1=0: V1=0 
FOR A1 = 0 TO 1.56 STEP .04 
A2 = ARCSIN(A*SIN(A1) ) 

Ml = A* COS(Al)/COS(A2) 

M2 = (l/A)*COS(A2)/COS(Al) 

B2 = 2*#PI*L*COS(A2) 

D1 = COS(B2)*COS(B2) 

D2 = SIN(B2)*SIN(B2) 

T3 = 1/(D1+((M1+M2)*(M1+M2)/4*D2)) 

H=INT(A1"286) : V= INT(T3*500) 

PLOT ( (H-Hl) , (V-Vl) ,HEX(FB) ) 

HI = H: VI = V 

NEXT A1 

PLOT (-H,-V,) 

NEXT A 
GOTO 210 

REM ; SUBROUTINE OPTION 3 

REM ; THIS SUBROUTINE PLOTS T3 VS A1 (CONST A), PARAM. L2 
INPUT "IMPEDANCE RATIO (R2/R1)=",A 
SELECT PRINT 213(64) 

PRINT " IMPEDANCE RATIO (R2/R1)=",A 

PRINT " T3 (0-1) VS AKO-89 DEG), PARAM L2O-0.5)" 

PLOT (10) (0,-50,"+") 

PLOT (10) (50 ,0 ,"+") 

PLOT (-500,0,) 

FOR L2 = 0 TO .5 STEP .1 

H=0: V=0: H1=0: V1=0 

FOR A1 = 0 TO 1.56 STEP .04 

LI = 1/ (2*SQR(1-A*A*SIN(A1)*SIN(A1) ) ) 
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960 A2 = ARCSIN(A*SIN(A1) ) 

970 Ml = A*COS(Al)/COS(A2 ) 

980 M2 = (l/A)*COS(A2)/COS(Al) 

990 B2 = 2*#PI*L2*L1*C0S(A2) 

1000 Dl= COS(B2)*COS(B2) 

1010 D2 = SIN(B2)*SIN(B2) 

1020 T3 = 1/ (Dl+( (M1+M2 )*(M1+M2 )/4*D2 ) ) 

1030 H = INT(A1*286): V = INT(T3*500) 

1040 PLOT ((H-Hl) , (V-V! ) ,HEX(FB) ) 

1050 HI = H: VI = V 

1060 NEXT A1 

1070 PLOT (-H,-V,) 

1080 NEXT L2 
1090 GOTO 210 

1100 REM ; SUBROUTINE OPTION 4 

1110 REM ; THIS SUBROUTINE PLOTS T3 VS L (CONST A), PARAM A1 
1120 INPUT "IMPEDANCE RATIO (R2/R10=", A 
1130 SELECT PRINT 213(64) 

1140 PRINT " IMPEDENCE RATIO ( R2 /R1 ) = " , A 

1150 PRINT " T3(0-l) VS L2(0-l), PARAMETER Al(22 DEG INC)" 

1160 SELECT PRINT 413 
1170 PLOT (10) (0,-50,"+") 

1180 PLOT (10) (50,0,"+") 

1190 PLOT (-500,0,) 

1200 FOR A1 = 0 TO 1.56 STEP .39 
1210 H=0: V=0: H1=0: V1=0 

1220 LI = 1/ (2*SQR(1-A*A*SIN(A1)*SIN(A1) ) ) 

1230 FOR L2 = 0 TO 1. STEP .02 
1240 A2 = ARCSIN(A*SIN(A1) ) 

1250 Ml = A*COS(Al)/COS(A2) 

1260 M2 = (l/A)*COS(A2 )/COS(Al) 

1270 B2 = 2*#PI*L2*L1*C0S(A2 ) 

1280 D1 = COS(B2 )*COS(B2) 

1290 D2 = SIN(B2 )*SIN(B2) 

1300 T3 = l/(Dl+((Ml + M2)*(Ml + M2)/4-’-D2)) 

1310 H = INT(L2*500): V = INT(T3*500) 

1320 PLOT ((H-H1),(V-V1) , HEX(FB)) 

1330 HI = H: VI = V 

1340 NEXT L2 

1350 PLOT (-H,-V,) 

1360 NEXT A1 
1370 GOTO 210 

1380 REM ; SUBROUTINE OPTION 5 

1390 REM ; THIS SUBROUTINE PLOTS T3 VS A (CONST L2 ) , PARAM A1 
1400 INPUT "BUBBLE LAYER/WAVELENGTH RATIO=",L 
1410 SELECT PRINT 213(64) 

1420 PRINT " BUBBLE LAYER/WAVELENGTH RATIO=",L 

1430 PRINT " T2(0-l) VS A(O-l), PARAM AK22 DEG INC) 

1440 PLOT (10) (0,-50,"+") 

1450 PLOT (10) (50,0,"+") 

1460 PLOT (-500,0,) 
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1470 

1480 

1490 

1500 
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1540 

1550 

1560 
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1720 
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FOR A1 = 0 TO 1.56 STEP . 39 
H=0: V=0: H1=0: V1=0 
FOR A = . 1 TO .9 STEP .04 
A2 = ARCSIN(A*SIN(A1) ) 

Ml = A*COS(Al)/COS(A2) 

M2 = (l/A)*COS(A2)/COS(Al) 

B2 = 2*#PI*L*COS(A2) 

D1 = COS(B2)*COS(B2) 

D2 = SIN(B2)*SIN(B2) 

T3 = l/(Dl+( (M1+M2)*(M1+M2)/4*D2) ) 

H = INT(A*500): V = INT(T3*500) 

PLOT ( (H-Hl) , (V-Vl) ,HEX(FB) ) 

H1=H: V1=V 
NEXT A 

PLOT (-H,-V,) 

NEXT A1 
GOTO 210 

REM ; SUBROUTINE OPTION 6 

SELECT R 

REM ; THIS SUBROUTINE PLOTS T3 VS A (CONST Al) , PARAM L2 
INPUT"ANGLE OF INCIDENCE= " , Al 
SELECT PRINT 213(64) 

PRINT "ANGLE OF INCIDENCE= " , Al 

PRINT "T3(0-D VS A(O-l), PARAMETER L2(0-.5, .1 STEPS)" 

PLOT (10) (0,-50,"+") 

PLOT (10) (50,0,"+") 

PLOT (-500,0,) 

FOR L2 = 0 TO .5 STEP .1 

H=0: H1=0: V=0: V1=0 

FOR A = .1 TO .9 STEP .04 

LI = 1/(2*SQR(1-A*A*SIN(A1)*SIN(A1) ) ) 

A2 = ARCSIN(A*SIN(A1) ) 

Ml = A*COS(Al)/COS(A2) 

M2 = (l/A)*COS(A2 )/COS(Al) 

L = L2*L1 

B2 = 2*#PI*L*COS(A2) 

D1 = COS(B2)*COS(B2) 

D2 = SIN(B2)*SIN(B2) 

T3 = 1/(D1+((M1+M2)*(M1+M2)/4*D2)) 

H = INT(A*500): V = INT(T3*500) 

PLOT ( (H-Hl) , (V-Vl) ,HEX(FB) ) 

HI = H: VI = V 
NEXT A 

PLOT (-H,-V,) 

NEXT L2 
GOTO 210 
END 
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